Hole states in wide band-gap diluted magnetic semiconductors and oxides 
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Puzzling disagreement between photoemission and optical findings in magnetically doped GaN 
and ZnO is explained within a generalized alloy theory. The obtained results show that the strong 
coupling between valence-band holes and localized spins gives rise to a mid-gap Zhang-Rice-like 
state, to a sign reversal of the apparent p-d exchange integral, and to an increase of the band-gap 
with the magnetic ion concentration. 
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Recent progress in the understanding of carrier- 
controlled ferromagnetic semiconductors such as 
Gai-a^Mn^^As and p-type Cdi-a^Mn^^Te [1] relies to a 
large extend on comprehensive optical, magnetooptical, 
photoemission, and x-ray absorption (XAS) studies, 
which provide mutually consistent information on 
exchange couplings between the band states and the 
localized spins as well as on the positions of d-bands 
with respect to the band edges. However, it becomes 
more and more apparent that this clear cut picture 
breaks down entirely in the case of nitride diluted 
magnetic semiconductors (DMS) and diluted magnetic 
oxides (DM0). Here, according to photoemission [I, 0] 
and XAS Q works, the relevant Mn^+/Mn^+ level 
resides in the valence band in both GaN [s*] and ZnO 
[2j, ^^], similarly to the case of arsenides, tellurides, and 
selenides. Furthermore, the evaluated magnitude of the 
p-d exchange energy PNq shows the expected enhance- 
ment (due to a stronger p-d hybridization [5] in nitrides 
and oxides) to the values of —1.6 eV in Gai-a^Mn^jN 
[S^ and to -2.7 eV |2|] or -3 eV [4] in Zni_^Mn^O, 
significantly over the value of —1 eV, typical for other 
DMS [1, 0, S]- contrast, optical and magnetooptical 
studies appear to reveal the presence of Mn^+/Mn^+ 
level in the band gap of both GaN [6] and ZnO [7]. 
Moreover, the determined values of /3No show either 
opposite sign and/or much reduced amplitude comparing 
to those steaming from photoemission and XAS as well 
as expected from the chemical trends. For instance, a 
detailed examination of the giant spin-splitting of free 
excitons in Gai-^^Mn^^N leads to PNq = +1.4 ± O.SeV 
[1], if the s-d exchange integral aNo = 0.2 ± 0.1 eV is 
assumed, though the probed charge state may differ from 
that examined by photoemission. Similarly, the study of 
magnetoluminescence of bound excitons in Zni-a^Mn^^O 
implies \PNo\ « 0.1 eV [9]. Interestingly, the contradic- 
tion in question is not limited to Mn-based nitrides and 
oxides - it is also rather evident for Zni_a^Coa;0, where 
PNq = —3.4 eV according to XAS and in agreement with 
the chemical trends while a thorough investigation 



of the free exciton splitting results in PNq ^ —0.6 eV or 
1 eV, depending on the assumed ordering of the valence 
band subbands [IQ] . 

In this paper we present arguments indicating that ni- 
trides and oxides belong to an original and unexplored 
class of DMS, in which the conditions for the strong cou- 
pling between the band hole and the localized magnetic 
impurity are met. We then calculate the spectral den- 
sity A{u;) of the band-edge states by a generalized alloy 
theory and show that in the strong coupling limit A{uj) 
exhibits two maxima corresponding to the bonding and 
anti-bonding states, respectively. We identify the former 
as the deep band-gap hole trap observed optically, while 
the latter as the position of the band-edge, and thus of 
the free exciton, renormalized in a non-trivial way by the 
hole coupling to the system of randomly distributed mag- 
netic impurities. As we demonstrate, our model makes 
it possible to explain the reversed sign and the reduced 
magnitude of the apparent /3No, as determined by the 
free exciton splitting, as well as the dependence of the en- 
ergy gap on the magnetic ion concentration x. In this way 
we reconcile the results of photoemission, XAS, optical, 
and magnetooptical studies and show their consistency 
with the chemical trends. Furthermore, our findings pro- 
vide a new piece of evidence for the inverse order of the 
valence subbands in ZnO as well as reconfirm that the 
local spin density approximation (LSDA) results in too 
high energetic positions of the localized d states. 

The description of exciton splittings in the works 
referred to above has been carried out within the 
time-honored virtual-crystal approximation (VGA) and 
molecular-field approximation (MFA), according to 
which the interaction of carriers and localized spins leads 
to bands' splitting proportional to the relevant bare ex- 
change integral and magnetization of the subsystem of 
magnetic ions. However, while this approach describes 
very well the giant exciton splittings in the tellurides 



ll\ it has already been called into question in the case 



of Gdi-a^Mn^^S [12]. In this system, unexpected depen- 
dencies of the band-gap and of the apparent exchange 
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integral (S^'^pp^ on X have been explained by circumvent- 
ing VGA and MFA either employing a non-perturbative 
Wigner-Seitz-type model [l2[ or by generalizing to DMS 
the alloy theory [1^. It has been found that the physics 
of the problem is governed by the ratio of a characteristic 
magnitude of the total magnetic impurity potential U to 
its critical value Uc at which a bound state starts to form. 
In particular, the weak coupling regime, where VGA and 
MFA apply, corresponds to U/Uc <^ I. By modeling 
the total potential of the isoelectronic magnetic impurity 
with a square well of the radius b one obtains [l2l.ll3|. 



U/Uc = 6m*[W - (5 + l)/3/2]/(7r''r6). 



(1) 



where the bare valence band offset WNq = dEy{x)/dx^ 
S is the impurity spin, and m* is the effective mass of 
a particle with the spin 5 = 1/2 assumed to reside in a 
simple parabolic band. 

In order to evaluate U/Uc for specific materials, we 
adopt p = —0.057 eVnm^ i.e., PNq ^ a^^ , as implied 
by the chemical trends [l, Furthermore, we note 

that the value of b lies presumably between the anion- 
cation and cation-cation distance, >/3ao/4 ^b < ao/V^. 
Hence, for Gdi_^Mn^Te, where S = 5/2, m*/mo = 0.65 
and WNo = -0.63 eV ^ we obtain 0.20 < U/Uc < 
0.33, in agreement with the notion that VGA and MFA 
can be applied to this system. In the case of Gdi-^^Mn^^S, 
where m*/mo = 1.0 and WNq = 0.5 eV pJi] the coupling 
strength increases to 0.77 < U/Uc ^ 1.25. Hence, in 
agreement with experimental findings, Gdi-^^Mn^^S rep- 
resents a marginal case, in which the local spin dependent 
potential introduced by the magnetic ion is too weak to 
bind the hole, but too strong to be described by VGA 
and MFA 

We argue here that even greater magnitudes of U/Uc 
can be expected for nitrides and oxides, as ao is smaller, 
while both m* and W tend to be larger comparing to 
compounds of heavier anions. We assume m* = 1.3mo 
and, in order to estimate W, we note that if the va- 



lence band offset is entirely determined by the p-d hy- 
bridization, W and P are related [13]. In the case of 
Gai-a^Mn^^N, for which the Hubbard correlation energy 
for the d electrons U^^^^^ = 10.4 eV and the position 
of the d-state with respect to the top of the valence 



band, e 



(e//) 



-5.7 eV [3] (determined so-far with accu- 



racy of about 1 eV), we arrive to 0.96 < U/Ucr ^ 1-6. 
Similarly for Zni-^Mn^^O, where U^'^ff^ = 9.2 eV and 
-1.5 eV [4], we find 2.0 < U/Uc < 3.3. In- 
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terestingly, the above evaluation of U/Uc for particular 
materials remains valid for DMS and DM0 containing 
other transition metals than Mn in the 2+ state as, to 
a first approximation, /3S does not depend on S [5]. Al- 
though the quoted values of U/Uc are subject of uncer- 
tainty stemming from the limited accuracy of the input 
parameters as well as from the approximate treatment of 
the hole band structure, we are in the position to con- 
clude that magnetically doped nitrides and oxides are in 
the strong coupling regime, U/Uc > 1, where the hole 
states cannot be described by VGA and MFA. 

We evaluate the effect of magnetic impurities upon a 
single band-edge particle {k = 0) with the spin s = 1/2 
by employing a generalized alloy theory developed by 
Tworzydlo [13]. The theory is built for non-interacting, 
randomly distributed, and fluctuating quantum spins S 
characterized by the field-induced averaged polarization 
{Sz{T, H)). The potential of the individual impurities is 
modeled [12] by the square- well potential containing both 
spin-dependent (exchange) and spin-independent (chem- 
ical shift) contributions, as discussed above. The par- 
ticle self-energy Ss^(cc;) is derived from the Matsubara 
formalism by summing up an infinite series for the irre- 
ducible self-energy in the average t-matrix approximation 
(ATA) [13]. Because of its general applicability for vari- 
ous alloys, we write down the form of (cj) explicitly, 
neglecting the direct effect of the magnetic field on the 
band states, 



S,, (co) = [(5 + 1 + 2s,{S,))^o{S) + (5 - 2s,{S,))^o{-S - 1)]/{2S + 1), (2) 

where So( J) = xNo{vU - miTm''UU' In'b^ /h^)- v = Airb^/S] U = -{W + jp/2)/v; U' = U - Eyc/{vNo); 

Eye = -xNo{s,{S,)f3^W); A = i{f^-f^')exp{if^')/{f^smf^'^m'cosf^y, k = b(2m'' E /h^f^^- k' = b[2m\E-U) /h^Y^^; 

I = {3(1 + A) + 3f^'\l + A) + 2if^'^ + 3[(i + /^O^ + ^f^'^ - 1)] cos{2f^') + 3i[A>:'(2i(l + A) + 1^') - 1] sm{2f^')}/{12f^'^). (3) 

I 

Here E = cj -\- Eye + i7, where 7 is an intrinsic life-time states, 

broadening of the k = state, taken here to be 3 meV. 1 1 

As H = Im (4) 

The quantity of interest is the spectral density of whose maxima provide the position of the band edge £^0 
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FIG. 1: Energy distribution (grey scale) of the spectral den- 
sity As^ {(jj) at the band-edge (k = 0) as a function of the 
strength of the coupling between the spin 1/2 particle and 
the system of unpolarized (a) and spin-polarized (b) mag- 
netic impurities. The magnitude of U/Uc is changed by the 
range of the impurities' potential b. 



in the presence of spin and chemical disorder as a function 
of {Sz{T,H)). Most of relevant experiments for magnet- 
ically doped nitrides and oxides has been carried out at 
low concentrations x < 3%, where effects of interactions 
among localized spins, neglected in the present approach, 
are not yet important. In this case, {Sz{T^H)) can be 
determined from magnetization measurements or from 
the the partition function of spin hamiltonian for the 
relevant magnetic ion at given temperature T and the 
magnetic field H. For higher x, it is tempting to take 
into account effects of the short-range antiferroniagnetic 
superexchange according to the standard recipe pj, 
i. e., via replacement of x by Xe// and T by T + Taf, 
where Taf{x^T) > and Xeff{x^T) < x. However, in 
the strong coupling case this procedure is rather inaccu- 
rate as aniferromagnetically aligned pairs contribute also 
to the band-edge shift. 

Figure 1(a) presents a grey scale plot of the spectral 
density As^{uj) as a function of the coupling strength 
in the absence of spin polarization, {Sz) = 0, so that 



the spin degeneracy is conserved, Ai/2(^) = ^-i/2(^)- 
In the weak coupling limit, U/Uc <C 1, the band-edge 
position £^0 is obviously given by —xNqW. For higher 
U/Uc a downward shift of Eq is visible, particularly 
rapid for U/Uc > 1, when the magnetic ion potential 
is strong enough to bind a hole, even if the magnetic 
ion is an isoelectronic impurity. Such an effect is also 
expected within the dynamic mean-field approximation 
[l6[. The small magnetic polaron formed in this way, 
reminiscent of the Zhang- Rice singlet, will be built also 
of k states away from the Brillouin zone center, as dis- 
cussed previously [I^- This band-gap hole trap should 
be visible in photoionization experiments d^ ^ d^ + h 
+ e, as indeed observed in both n-Gai-a^Mua^N [6.] and 
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FIG. 2: Dependence of the band-edge position £o(a) and of its 
spin splitting As (b) on the concentration of magnetic impu- 
rities X for three values of the coupling strength in comparison 
to the values expected within VGA and MFA. 



Zni_a^Mna;0 ^] but assigned in those works to d^ ^ d^ + 
e transitions. In Gai-a^Mn^^N samples, in which donors 
were partly compensated, intra-center excitations corre- 
sponding to this level were also detected at 1.4 eV 
and analyzed in considerable details [H, [2^. While 
symmetry considerations cannot discriminate between 
d^ + h and d^ many-electron configurations, the larger 
crystal-field splitting and the smaller Huang-Rhys factor 
in Gai_a;Mn^N comparing to their values in (II,Cr)VI 
compounds 1^, 2^ point to a relatively large localization 
radius, expected for the d^ + h case. Our interpretation 
reconfirms also the limited accuracy of ab initio compu- 
tations within LSDA in the case of transition metal insu- 
lators, which place the d"^ Mn level within the band-gap 
of GaN [Hj. 



Interestingly, as seen in Fig. 1(a), when U/Uc in- 
creases beyond 1, the spectral density is gradually trans- 
ferred from the bonding state discussed above to an anti- 
bonding state appearing above the band-edge energy Eq 
expected within VGA and MFA. We identify this state 
with the actual valence band-edge for U/Uc > 1, whose 
position £^0 determines, e. g., the onset of interband op- 
tical transitions and the free exciton energy. In Fig. 2(a), 
we plot Eo{x) for U/Uc < I and U/Uc > 1 in compar- 
ison to the VGA and MFA expectations. We see that 
if U/Uc < 1, the correction to VGA and MFA leads 
to a reduction of the band-gap with X ^ as observed in 
Gdi_^Mn^S [13], Zni_^Mn^Se [22], and Gai_^Mn^As 
(1^. In contrast, for U/Uc > I the present model predicts 
an increase of the gap with x, in agreement with the data 
for Gai_^Mn^N [20], Zni_^Mn^O [23^, and Zni_^Go^O 

As shown in Fig. 1(b), the saturating magnetic field, 
such that {SziT^H)) = —2.5, produces a downward 
and upward shift of the bonding and anti-bonding states 
(which correspond to Sz = 1/2), respectively, and leads 
to the appearance of a non-bonding state {sz = —1/2), 
whose energy is virtually independent of U/Uc Re- 
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FIG. 3: Dependence of the band-edge spin splitting on the 
spin polarization, As{{Sz)) for three values of the coupling 
strength in comparison to the values expected within VGA 
and MFA, Aj'^^((5'^))(a). Normalized apparent p-d ex- 
change integral defined as As(— 2.5)/AJ^^^(— 2.5) for various 
combinations of the materials parameters (b) 



markably, the latter resides below the anti-bonding level, 
which means that the character of the apparent p-d ex- 
change changes from antiferromagnetic for U/Uc < 1 to 
ferromagnetic for U/Uc > 1. Extensive studies of the gi- 
ant free excitons' splitting Ag in Gai-a^Mua^N [8] and 
Zni_a^Coa^O [10] demonstrated a linear dependence of 
As on X and {Sz{T^H))^ as expected within VGA and 
MFA. Figures 2(b) and 3(a) demonstrate that approxi- 
mately linear dependence of As on x and {Sz{T^H)) is 
predicted within the present theory, too, except for an 
anticrossing behavior occurring when the non-bonding 
state is in a resonance with the non-normalized band- 
edge, ^0 = 0. Accordingly, for the sake of comparison 
to the experimental determinations. As can be charac- 
terized by an apparent p-d exchange integral according 
to A, = -xNoP^^'pp^SzIt.H)). In Fig. 3(b), we depict 
the expected evolution of (3^^^^^ / (3 with U/Uc for several 
values of W and m* . We see that our theory explains the 
sign reversal of [3^^^^^ and its reduced amplitude compar- 
ing to the P values determined by the photoemission and 
XAS experiments. Our findings may imply also that the 
exchange coupling between the band hole and the hole 
on the d^ + h center is not qualitatively important in 
Gai-a^Mn^^N as well as verify the reversed order of va- 
lence band subbands in ZnO [24]. 

In summary, the approach put forward here allows one 
to reconcile the findings of photoemission, XAS, optical, 
and magnetooptical measurements accumulated over re- 
cent years for the magnetically doped nitrides and oxides. 
In view of our results, these systems form an outstanding 
class of materials, in which a number of concepts devel- 
oped earlier for DMS has to be revised. In particular, 
the strong hole localization discussed here, along with 
the issues of self-compensation and solubility, need to be 
overcame in order to observe the hole-mediated ferromag- 
netism in nitrides and oxides. 
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